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We consider a thin superconducting film with randomly magnetized dots on top of it. The 
dots produce a disordered pinning potential for vortices in the film. We show that for dots with 
permanent and random magnetization normal or parallel to the film surface, our system is an 
experimental realization of the two-dimensional XY model with random phase shifts. The low- 
temperature superconducting phase, that exists without magnetic dots, survives in the presence of 
magnetic dots for sufficiently small disorder. 
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I. INTRODUCTION 

Since the pioneering papers by Berezinskii 1 and 
Kosterlitz and Thouless^ finite temperature phase tran- 
sitions in two dimensional (2D) systems, which have a 
continuous symmetry specified by a phase, have become 
a very active research field. It turned out that two di- 
mensional superfluids 3 - as well as thin superconducting 
films^— have a Berezinskii-Kosterlitz-Thouless(BKT) 
transition at a finite temperature Tib- The low- 
temperature phase is superfluid (superconducting) and 
has quasi- long range order, while the disordered high- 
temperature is normal liquid (metallic) for superfluids 
(superconductors). These systems are successfully de- 
scribed by the 2D XY model, see Eq. (P with Ay = 0. 
Theoretical predictions were confirmed experimentally: a 
universal jump in the superfluid density at Tirr and dif- 
ferent current- voltage characteristics of superconducting 
films below and above T^d^- 

Introducing disorder through random phase shifts in 
the 2D XY model physics becomes more complex. The 
Hamiltonian of the 2D XY model with random phase 
shifts readsiS^- 
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where the sum runs over all nearest neighboring sites on a 
square lattice. <f>i denotes the phase of the order parame- 
ter, while Ay are quenched random phase shifts on bonds 
of the square lattice produced by some kind of disorder. 

We ctSSUH16 t licit -^ij 

are Gaussian distributed, uncorre- 



cted, have the zero mean (Aj 



and the variance 



(Afj) = a. The pure 2D XY model contains two different 
kinds of excitations which are decoupled. The spin-wave 
excitations describe small changes of the phase <j> and do 
not drive any phase transitions. Another kind of excita- 
tions are vortices, topological excitations, and they are 
essential for the existence of a BKT transition in the 2D 
XY model at T2D- Vortices are bound in pairs in the low- 
temperature phase, and unbound at high temperatures. 
Introducing the disorder in the form of random phase 
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FIG. 1: Phase diagram of the 2D XY model with ran- 
dom phase shifts. At low temeprature and weak disorder 
a < o c a superconducting phase (SC) occurs, and there vor- 
tices are bound in pairs. Above T2D vortices are unbound 
and the phase is non-superconducting (N). According to some 
studiesi£~-i2 at low temperatures and for strong enough disor- 
der appears another superconducting phase. 



shifts, the low-temperature superconducting phase sur- 
vives for weak enough disorder, as has been shown first 
analytically^— and then numericall y 15 ' 16 . The phase 
diagram of the model ((TJ is given in Fig. Q] 

The Hamiltonian ([T]) describes the thermodynamic 
behavior of several disordered systems, including two- 
dimensional ferromagnets with random Dzyaloshinskii- 
Moriya interaction^, Josephson junction arrays with po- 
sitional disorder— and vortex glasses 2 ^. In this paper we 
will consider another system and show that it also be- 
longs to the class of systems that can be described by 
the model (JTJ). It is a ferromagnet-superconductor hy- 
brid system. The system consists of a thin supercon- 
ducting film covered by magnetic dots with permanent, 
but random magnetization. For more details about hy- 
brid systems see recent review s 22 ' 23 . We will show that 
our system can be described by the Hamiltonian (JTJ un- 
der some conditions defined below. Knowing the solution 
of the model ((T|) we can establish possible phases of the 
system. 

Placed on top of the film, a single magnetic dot with 
sufficiently large magnetization, normal to the film sur- 
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face, induces and pins vortices and antivortices around 
i^24-^27_ A regular lattice of magnetic dots with con- 
stant and equal magnetization is the source of a peri- 
odic pinning potential for vortices in the film. The mag- 
netic dots induce periodic arrays of vortices and antivor- 
tices in the film, provided the dot magnetization is high 
enoug h 28 ! 29 . The periodic pinning is absent in the case 
of a random dot magnetization, and there magnetic dots 
are a source of disordered pinning potential for vortices. 
Lyuksyutov and Pokrovsky 22 i 24 considered a supercon- 
ducting film with a magnetic dots array with random, 
sufficiently strong magnetic moments and concluded that 
the dot array induces the resistive state in the film. They 
have not considered the case when magnetic moments 
are weak. In this paper we will show that the supercon- 
ducting state survives provided the disorder is not too 
strong. Hybrid systems with a regular lattice of mag- 
netic dots with random magnetization 3 ^, and without a 
lattice, but homogeneously distributed dots on top of a 
thin superconducting film^i have been recently experi- 
mentally realized. 

In the rest of the paper in section II we introduce the 
model for our hybrid system and give its solution by map- 
ping it to the model JT]). Section III contains discussions 
and conclusions. Some technical details are postponed to 
the appendix. 



II. MODEL AND ITS SOLUTION 

We consider a thin superconducting film character- 
ized by the London penetration depth A^, the coherence 
length £, a thickness d and a typical lateral dimension L. 
The London penetration depth is in general temperature 
dependent and we assume that our film has the effective 
penetration depth A = A 2 /d that exceeds film's lateral 
dimension L. This limit is valid for "dirty" supercon- 
ducting films 4 . In addition, provided the bulk critical 
temperature is larger than a critical temperature T 2 d 
for vortex unbinding, the film has a BKT transition at 

We consider the dots with permanent random magneti- 
zation placed on top of the film. They produce a random 
potential V for vortices in the film. Assuming vortices of 
vorticities rij are places on a quadratic lattice with the 
lattice constant a (which is of the order of the coherence 
length), the effective lattice Hamiltonian for the system 
may be written as 

H v = ^ [ n K E c + U v ) + riiVi] + \ y^ j n l n j U vv (p ij ) 1 



(2) 



where the sum runs over all lattice sites, £0 = 
0q/(I67t 2 A), the flux quantum is 0o = hc/(2e), Vi is the 
random potential at site i, while E c is the single vortex 
core energy which is of the order e - U v and U vv (pij) 
are the single vortex energy and the interaction energy 



of two vortices separated by a distance p^ respectively, 
and for L < A read 3 ^ 3 . 

U v = e In — , U vv (pij ) = 2e In — . (3) 

Q Pij 

Using the expressions <j3j) it is useful to rewrite the Hamil- 
tonian ^ in the form 



i i^j 
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where we have introduced the total vorticity of the sys- 
tem N — X)i n 2- I n the limit L 3> a and without the 
disorder potential Vi — 0, the last term in Eq. ((4]) penal- 
izes the total energy for nonzero total vorticitie o 34 ' 35 , so 
one has N = 0. Then, a superconducting film described 
by the model ([4]) has a BKT transition at the tempera- 
ture T 2 d = £o(?2d)/2, where £0(^20) denoted that one 
should take value X(T2d) renormalized by the presence 
of vortices2^£. 

Next, we consider effects of the random potential. First 
we consider the case when the dots have random mag- 
netization parallel to the film surface. We model them 
as magnetic dipoles placed on top of the film at lattice 
sites. To characterize statistical properties of the dots, 
we assume that the x— and y— components of the mag- 
netic moment at site i are Gaussian distributed, have 
zero mean value and are uncorrelated from site to site: 

(mia) = 0, (mi a mjp) = M 2 5ij5 a p, a, (3 — x,y (5) 

where (. . .) denotes an average over disorder. Since the 
dots have random magnetic moments, Ai is the measure 
for a typical magnetic moment of a magnetic dot at some 
lattice site. 

The interaction energy between a single dot having the 
magnetic moment m parallel to the film surface and a 
vortex of vorticity n placed at a relative distance p from 
the dot can be calculated using the approach developed 
in Ref. 36 , and reads 2 ^ 



Ul v (n,p) 



m • p n<po 
p ~2tt 

m • p ncj)Q 
p I6A 2 



dk 



kJi(kp) 



1 + 2Xk 



(6) 



where J\ is the Bessel function of the first kind, Hi is the 
Struve function and Y\ is the Bessel function of the sec- 
ond kind 37 . Notice that the interaction energy between 
the magnetic dipole and the vortex can be simply writ- 
ten in the form — m • B„, where Y5 V is the magnetic field 
produced by the vortex at the dipole positio n 27 ' 38 . In the 
limit p <C A the interaction energy (jJJ]) reads 



UL v (n,p) = 



m • p n0o 

p AirXp 



(7) 
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The random site potential Vi is given as a sum over the 
lattice 



Vi 



4ttA ^ 



Pi 



(8) 



By summing over j : ^ i in the previous expression we 
avoid the short scale cutoff divergence of the interaction 
energy ([6]) at p = 0, which exists because the dots are 
placed at the top of the film in our model. In reality 
magnetic dipoles are separated from the film surface by 
some small distance ~ a. Since Vi is a sum of many 
independent random variables it is Gaussian distributed 
(notice here that the assumption about Gaussian distri- 
bution for mi a is not necessary condition for V% to be 
Gaussian distributed; it is sufficient that dots have a dis- 
tribution with a finite variance). Its mean, variance and 
site to site correlations can be calculated and read 
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The model Q with the disorder potential (JSJ which 
has properties (HJ, (fT0|) and ([TT]) matches the vortex part 
of the 2D XY model with random phase shifts—^. From 
the solution of the model (TT J 11 ! 14 ! 39 we know its phase di- 
agram, see Fig. [1] At zero temperature there is a critical 
value for the typical magnetic moment per unit length 



Mc = 
a 8irV2n 



(12) 



For Ai < Ai c the system has no free vortices N = 0, 
while for M. > A4 C the disorder spontaneously creates 
and induces unbound vortices, and one generally has 
W / 0. A simple argument — based on Eq. Q which 
compares the energy loss for the single vortex creation 
and the energy gain due to disorder fluctuations also gives 
the critical value (|12[) for zero temperature. The connec- 
tion between a introduced as the disorder strength in 
Eq. ([T]) and the typical magnetic moment A4 introduced 

in Eq. ([5} is a = ( 4 ^ ^ ^ , while the random phase is 



A, 



4tH 
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e z x mi where the x— (y—) component of Aj 



corresponds to the disorder Aij on horizontal (vertical) 
bond at site i. 

In the case that the dot magnetization is normal to 
the film surface, similar to Eq. ([6]), for the interaction 
energy between a dot of magnetization m and a vortex 
of vorticity n separated by a distance p we get 



U^ v {n,p) = m 
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The previous expression simplifies for p <C A and reads 

n4> 



U mv (n,p) = m 



AnXp' 



(14) 



By assuming that the film is covered by magnetic dots 
with random magnetization normal to the film surface 
(along the z direction), satisfying 



(m iz ) = 0, (m lz m jz ) = M 2 6ij, 
for the site random potential we get 

i ; ''j 4ttA^ pij 



(15) 



(16) 



which, as we show in the appendix, is equivalent to Vi 
and hence has statistical properties (|9|)- (fTT|) . We may 
conclude that the system with magnetic moments normal 
to the film has the same phase diagram as in the case of 
moments parallel to the film surface. 



III. DISCUSSIONS AND CONCLUSIONS 

Having shown the equivalence between our system and 
the vortex part of the 2D XY model with random phase 
shifts, we may infer some properties of the former. The 
phase diagram is given in Fig. [TJ The low temperature 
and low disorder phase is superconducting. There vor- 
tices and antivortices are bound in pairs. The current- 
voltage characteristic for T — > T2D and T < T2D and for 
weak disorder is expected to be very similar to the one of 
Halperin and Nelsoni for the pure case, V ~ I 3 , possibly 
with a small correction due to the disorder. This phase 
has zero linear resistivity. The high temperature phase is 
metallic and has a nonzero linear resistivity. There free 
vortices dissipate energy and produce the linear current- 
voltage characteristic V ~ I. 

In Refi^ the conclusion about the resistive state of 
the film when the dots with normal magnetization are 
present relies on the assumption that the randomly mag- 
netized dots pin vortices of vorticity ±1. These pinned 
vortices serve as a source of the random potential for 
other bound vortices, which unbind and fill deep valleys 
of the random potential. These unbound vortices lead 
to the resistive state of the film. We agree that this 
scenario occurs for sufficiently strong disorder when the 
dots can induce and pin vortices. A single dot can induce 
and pin quite different configurations of vortices and an- 
tivortices regarding its magnetic moment- 7 . We expect 
that a random lattice of dots can also pin, from site to 
site, quite a different number of vortices and antivortices 
which produce different potential than one assumed in 
Refi^. However, our conclusion, that the resistive state 
in films occurs when the disorder is sufficiently strong, 
agrees with the one from Ref^. Moreover, we give the 
strength of the disorder above which the resistive state 
occurs. 
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By making a comparison between A4 C and the mag- 
netic moment m lc of a single magnetic dot with normal 
magnetization necessary to induce and pin an extra vor- 
tex in the film, we obtain m\ c rs M c \f&K{b/a) ln(L/a), 
where b is the distance between the dipole and the film 
surface. In addition, knowing that the value M c corre- 
sponds not to typical but rare magnetic moments from 
the tail of distributio n 11 ! 40 , we may conclude that even a 
very rare magnetic dot in the film that has the magnetic 
moment A4 C is not able to induce and pin vortices. Such 
pinned vortices served as a source of random potential in 
Ref^i. 

In this paper we have considered the dots as magnetic 
dipoles. This fact is unimportant as long as the dot size is 
not too big with respect to the lattice constant. What is 
crucial for any kind of magnetic dots is their interaction 
with a vortex which decays as 1/p, which is universal for 
any geometrical shape of dots, when the vortex is suffi- 
ciently far from the dot. The magnetic field produced by 
a vortex decays as 1/p and the interaction energy dot- 
vortex universally decays, regardless of the shape of the 
dot. This form of the interaction produces logarithmi- 
cally diverging, with the system size, variance of the dis- 
order potential that is characteristic for the Hamiltonian 

Recently, the question of a possible third phase for 
strong disorder and at low temperatures has been raised 
in numerical studies of the model ([1]) with uniformly ran- 
dom phase Aij in the interval [— r7r, r7r] with < r < 1, 
with conflicting results about its existenc o 17 ' 19 ! 41 . While 
for r < r c tx 0.37 16 it is accepted that the superconduct- 
ing phase survives at low temperature, the case r > r c 
is still under debate. An experimental study of Joseph- 
son junction arrays with positional disorder 18 supports 
the existence the third phase. The new phase is super- 
conducting according to the experiments of Ref. 1 — and 
numerical investigations of Ref. 1 —. In the limit a — > oo 
the model ((TJ and the so-called gauge glass limit, r = 1, 
are equivalent, and we expect that the conjectures about 
the phase diagram for the case of uniformly distributed 
phase apply as well to the Gaussian distribution of the 
phase. The existence of the third phase is a challenging 
question that could be experimentally resolved by using 
magneto-superconducting hybrid systems as a realization 
for the model ([1}. 

To conclude, we have shown that a thin superconduct- 
ing film covered by magnetic dots with random magneti- 
zation provides an experimental realization for the two- 
dimensional XY model with random phase shifts. The 
phase diagram of the latter model helped us to conclude 
that a low-temperature superconducting phase of a su- 
perconducting film without dots survives when the dots 
are placed on top of the film, provided their magnetiza- 
tion is not too large. 
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IV. APPENDIX 

In this appendix we prove that expressions for the dis- 
order potential produced by magnetic moments paralel 
to the film surface ([5]) and normal to the film surface 
(Unj) are equivalent. Rewriting the expression va.j ■ p^/ p 1 - 
from Eq. (J8)) as rrij cos(o>j — a p )/pij, where coso, = 
m.j ■ e x /rrij and cosa p = p^ ■ e x /pij, we will in the 
following show that the distribution of random variable 
m r j = rrij cos(<x/ — a p ) (which is the projection of nx, 
onto is Gaussian, with zero mean and the variance 
M 2 . 

By assumption ([5]), the components of the magnetic 
moment m; x and Gaussian distributed and have 

the distribution function 

P(t) ^V^M eXP (^) t = m *> m »- (17) 
Then the distribution function of the magnetic moment 

™3 = \f m % + m2 jy 1S 

/OO pOO 
-oo J — oo 

x6 (m 3 - - yfinf x + m 2 y ) = ^ exp (~^J , (18) 

while the angle ccj between m,j X and rrij is uniformly dis- 
tributed in the interval [0, 2n) and has the distribution 
p{ocj) = l/(27r). The distribution of the random variable 
m r j is 

p(m r j) = / drrij / dajp(m.j)p(aj) 
Jo Jo 

xS (m r j — rrij cos(a 3 — a p )) . (19) 

The previous expression can be most easily evaluated first 
by taking the Fourier transform of p(m r j) 

Pk = J dm iJ da ^W eXP (-w) 

x exp[ikmj cos(aj — a p )] = exp ( — J , (20) 

and then taking the inverse Fourier transform of the pre- 
vious expression. It leads to 

1 / m 2 \ 

Km - )= V^ eXP V2^J- (21) 
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In that way we have proved that the random potential 

Vi -^\2^ pi - ^ x 2^ p w 



matches the random potential (|T6|) . We conclude that 
frozen magnetic dipoles parallel to the film create the 
same random potential for vortices in the film as mag- 
netic dipoles normal to the film, provided both are Gaus- 
sian distributed. 
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